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BUCKLING  OF  CLAMPED  SHALLOW  SPHERICAL  SHELL: 


by 

Bernard  Bud Lansky 


L.  INTRODUCTION 

This  paper  is  concerned  with  the  axisymmetric  buckling  under  uniform 
pressure  of  a  shallow  portion  of  a  spherical  shell,  clamped  along  a 
circular  boundary  (Figo  1).  Toe  same  problem  has  been  the  subject  of 
several  previous  theoretical  investigations  (Refs.  l-7)«  Although  these 
previous  studies  were  all  based  on  the  same  non-linear  differential 
equations,  the  results  found,  by  a  variety  of  techniques,  are  generally  in 
disagreement  with  each  other,,  Some  of  the  trouble  in  these  solutions  is 
evidently  related  to  the  waviness  of  the  shell  distortions,  which  tends  to 
increase  with  decreasing  shell  thickness.  The  present  analysis,  winch 
exploits  an  integral -equation  formulation  of  the  problem,  should  give 
results  of  uniform  accuracy  for  all  thicknesses. 

Same  qualitative  information  concerning  the  buckling;  of  complete 
spherical  shells  may  be  recalled.  According  to  the  classical  theory,  a 
complete,  initially  perfect  spherical  shell  subjected  to  external  pressure 
contracts  in  size,  without  bending,  until  the  buckling  pressure  q^  is 
reached;  after  buckling  (Ref.  8),  the  pressure  needed  for  equilibrium 
drops  sharply  (see  curve  A,  Fig.  2).  On  the  other  hand,  initial 
geometrical  imperfections  in  the  complete  spherical  shell  permit  bending 
to  start  'mmediately  upon  the  application  of  pressure,  and  the  subsequent 
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non-lineor  variation  of  pressure  with  deflection  leads  to  a  local  maximum 
of  the  pressure  (see  curves  B,  G,  Fig.  ?) ,  at  which  point  buckling  - 
or  snapping  -  occurs.  Experiments  show  that  buckling  pressures  of  complete 
spherical  shells  are  substantially  lower  than  ,  arid  it  is  now  generally 
believed  that  initial  imperfections  are  responsible  for  this  disagreement. 

In  contrast  to  the  behavior  of  the  perfect  complete  sphere,  the 
perfect  clamped  shallow  shell  distorts  axisymme trie ally  out  of  the 
spherical  shape  at  the  beginning  of  loading;  in  effect,  the  pressure  makes 
its  owT'  "  initial"  imperfections,  and  so,  as  in  the  case  of  the  imperfect 
complete  sphere,  buckling  of  the  perfect  clamped  shell  may  be  related  to 
the  local  maximum  of  the  pressure  an  bending  proceeds.  It  does  not 

necessarily  follow,  however,  that  initial  geometrical  imperfections  in  the 
shallow  shell  may  be  disregarded.  In  fact,  for  sufficiently  small  shell 
thickness,  the  early  bending  deformation  of  a  perfect  shallow  spherical 
shell  under  uniform  pressure  is  confined  to  a  narrow  "  boundary  layer" 
near  the  edge,  with  the  rest  of  the  shell  remaining  essentially  spherical 
(Ref.  9).  Consequently,  as  thickness  decreases,  q^  of  the  initially 
perfect  shallow  shell  is  expected  to  approach  q^  of  the  corresponding 
complete  sphere,  and  experimental  confirmation  would  then  not  be 
anticipated  unless  geometrical  .^perfections  were  taken  jnto  account. 

In  the  present  paper,  results  are  found  for  the  axisymmetric  buckling 
pressures  of  perfect  clamped  shells,  as  well  as  of  shells  having  small 
axis ymme trie  initial  geometric  imperfections.  (Similar  imperfections 
were  also  considered  in  Ref.  7).  It  should  be  noted,  however,  that 
axisymmetric  buckling  does  not  necessarily  correspond  to  the  lowest 
buckling  pressure  of  the  initially  perfect  shell;  the  possibility  exists 
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that  an  instability  associated  with  unsymmetrical  distortions  may  occur 
before  the  axisymmetric  buckling  pressure  is  reached.  Furthermore,  un¬ 
symmetrical  initial  imperfections  might  have  a  more  significant  effect 
than  axisymmetric  ones;  although  this  conjecture  is  not  explored 
analytically  in  this  paper,  some  qualitative  observations  concerning  its 
likelihood  are  made. 


2.  BASIC  EQUATIONS 

The  basic  relations  of  non-linear  shallow  shell  theory  (Ref.  10) 
will  be  reproduced  for  the  case  of  axisymmetric  deformations  of  shells 
of  revolution  subjected  to  uniform  pressure. 

Moment  equilibrium  (see  Figs.  1  and  3)  requires 

4-(rM  )  -  M  -  rQ  =  0  (3) 

ar  r'  t 

horizontal  force  equilibrium  gives 


'  -  Nt 


and  vertical  equilibrium  of  a  polar  cap  gives 


(2) 


where  M  ,  M.  ,  N  ,  N, 
r  ’  t  r  t 

length  as  shown  in  Fig.  3 
ward  displacement,  and  z 


rQ  SSL.  +  rN  _  -~)  =  0 
2  r  dr  dr 


dz 


(3) 


and  Q  are  moments  and  forces  per  unit 
,  q  is  the  pressure,  w  is  the  vertical  down- 
defines  the  shell  middle  surface,  as  shown 


in  Fig.  1. 


The  generalized  force -distortion  relations  are 
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where  u  is  the  horizontal  displacement,  t  is  the  shell  thickness. 

E  is  Young’s  modulus,  and  v  is  Poisson's  ratio. 

Helating  the  membrane  forces  aid  to  a  stress  function  y 

by 


N 


and 


=  1 


dr 


(8) 


(?) 


provides  satisfaction  of  (2);  eliminating  M  ,  M,  and  Q  from  (1), 


(3),  (6),  and  (7)  gives 


EtJ 


12  (1-v  )  | 


da 


s(r "  SI  ~  (^}t  =  _qr  +  w  (10) 


dr 


dw 


where  p  =  -  —  is  the  rotation  of  an  element  of  the  shell  (see  Fig.  3). 
Finally,  elimination  of  the  displacement  u  from  (U)  and  (5)  provides 
the  compatibility  equation 

dz 


it  \pir  i£}  ■  r  j  +  (dT)p  =  -  \ p2  (11) 


The  simultaneous  equations  (10)  and  (11)  for  the  variables  p  and 
\|r  are  equivalent  to  the  equations  governing  the  analyses  of  Kefs.  1-7. 
The  boundary  conditions  appropriate  for  the  clamped  shall  are 


P  =  0  at  r  =  a 


(12) 
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aixd,  as  a  consequence  of  Eq.  (5)  and  the  requirement  that  u  vanish  at 
the  edge, 

8  "  V  I  =  0  at  r  "  a  (13) 

In  addition,  conditions  of  finite  ness  of  (3  and  f  are  imposed  at 

r  =  0  . 


Within  the  framework  of  shallow  shell  theory,  an  initial  spherical 
shape  is  represented  by  the  parabola 


z  =  H 
o 


(111) 


where  (see  Fig.  1)  H  is  the  shell  11  rise*  .  Then,  introduction  of  the 
geometrical  parameter  X  defined  by 

\h  =  U8(l-v2)(|)2  (15) 

and  substitution  of  the  non-dimensional  variables 


into  Eqs.  (10)  and  (11)  give  for  the  initially  perfect 
the  non-dimensional  equations 

(x3')'  -  ~  +  =  -2px  +  9$ 

(*$')'  -  |  -  x9  =  -  |  e2 


spherical  shell 

(19) 

(20) 


where  primes  denote  differentiation  with  respect  to  x  .  Here  the  non- 
dimensional  pressure  parameter  p  is  defined  as 


(21) 


where 
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P 


=  _a. 


q 


o 


(i)2 


(22) 


V3(l-v2) 

Is  the  classical  buckling  pressure  of  a  complete  spherical,  shell  having 


the  same  radius  of  curvature  R  =  —  as  the  shallow  shell  under  con- 

£ii 

sideration.  The  boundary  conditions  (12)  and  (13)  transform  to 


and 


e(x)  -  o 

(23) 

x$  »(x)  -  v#(X)  =  0 

(2U) 

Note  tijat  in  the  present  non-dimensional  formulation  the  geometrical 
parameter  \  determines  the  interval  in  the  x  coordinate  in  which  the 
differential  equations  (19)  and  (20)  apply,  but  does  not  appear  explicitly 
in  these  equations. 

According  to  Reissner  (Ref.  9)  the  shallowness  requirement  is  roughly 
H/a  =  1/6  .  This  restriction  is  quite  independent  of  the  parameter  X  , 
which  may  be  as  large  as  desired.  To  avoid  confusion,  it  is  best  to  regard 
H  and  a  as  fixed  in  magnitude,  and  to  think  of  small  X’s  as 
representing  thj.ck  shells,  whereas  large  X's  correspond  to  thin  shells. 

The  remark  in  Section  1  concerning  the  expected  variation  of  buckling 
pressure  with  thickness  can  now  be  made  more  precise:  as  X  becomes 
infinite,  Pcr(=  o,, )  for  the  initially  perfect  shallow  shell  should 
approach  1  . 

The  surface  of  a  shallow  shell  whose  initial  configuration  deviates 
axis ymme trie ally  from  perfect  sphericity  may  be  represented  by 
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*0  -  H  [l  -  (f)2  “  ee(r)]  (25) 

where  e(r)  is  the  shape  of  the  imperfection,  In  this  case,  the  ecyiations 
analogous  to  (19)  and  (20)  become 


e(r)  =  [l  -  (J)2  2  (29) 

For  this  choice,  e  equals  the  ratio  of  the  downward  initial  displace¬ 
ment  at  the  center  of  the  shell  to  the  rise  H  of  the  perfect  shell. 

It  should  be  noted  that,  in  shallow  shell  theory,  an  initial 
i arabolic  shape  does  not  constitute  a  deviation  from  "  perfect " 
sphericity;  thus  the  non-linear  differential  equation.:  (19)  and  (20) 
are  satisfied  exactly  by  the  membrane  state 

fr° 

[f  =  -2px 

This  state  does  not,  however,  satisfy  the  fixity  condition  (2U),  and  so 
bending  of  the  clamped  shell  must  oo.cur  for  all  p  . 
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3.  REVIEW  OF  PREVIOUS  WORK 

Figure  U  summarizes  the  results  of  previous  attempts  to  calculate 
Pcr  of  the  initially  perfect  clamped  shell*  Feodosiev  (Ref.  1.)  applied 
the  Galerkin  method,  using  simple  polynomial  functions,  to  obtain 
approximate  solutions  of  the  governing  non-linear  equations*  Kaplan  and 
Fung  (Ref.  2)  obtained  the  first  two  terns  in  a  power  series  expansion 
of  p  versus  w(0)  j  in  their  procedure,  the  non-linear  differential 
equations  were  reduced  to  a  set  of  linear  equations  by  a  perturbation 
technique.  They  were  able  to  find  results  for  ,  as  determined  by 

the  condition  =  0  s  only  for  relatively  low  values  of  X  . 

(Kaplan  and  Fung  also  obtained  one  improved  value  of  p^^  for  X  =  U 
by  solving  several  more  of  their  perturbation  equations  —  by  finite 
difference  approximations  —  to  get  four  terms  in  the  power  series  for 
p  )*  Archer  (Ref*  3)  recognized  that  for  higher  values  of  X  the 
variation  of  w(0)  with  p  is  not  necessarily  monotonic,  and 
accordingly  used  the  deflection  at  some  intermediate  station,  when 
necessary,  in  a  perturbation  expansion  similar  to  that  of  Kaplan  and 
Fung.  Archer  solved  his  linear  perturbation  equations  approximately 
by  means  of  a  finite  number  of  terms  in  a  Fourier-Bessel  expansion. 

Note  that  Archer's  results  for  p^_  do  appear  to  approach  1  as  X 
increases.  Reiss,  Greenberg,  and  Keller  (Ref*  h)  tried  power  series 
solutions  of  the  non-linear  differential  equations $  for  a  given  value 
of  X  ,  they  looked  for  solutions  at  successive  increments  of  p  ,  and 
presumed  that  p  had  been  exceeded  when  the  search  proved  fruitless. 
Their  results  are  embodied  in  the  collection  of  vertical  intervals 
shown  in  Fig.  k,  within  each  of  which  pc^  is  supposed  to  lie.  A  very 


-9- 


sinrLlar  power  series  technique  was  utilized  independently  by  Weinitschke 
(Ref,  5)  who  obtained  definite  results  for  pcr  only  for  X  -  5  >  but  at 
several  higher  dues  •  ‘  (see  the  dotted  curve  in  Fig.  <4),  he  succeeded 
in  finding  solutions  to  the  non-linear  equations  at  values  of  p  above 
those  where  Reiss,  Greenberg  and  Keller  failed.  According  to  Vfeinitschke, 
the  values  of  p  along  the  dotted  portion  of  his  curve  constitute  only 
lower  bounds  to  the  correct  values  of  pcr  .  Finally,  Chen  (Ref.  6), 
improving  upon  a  procedure  used  by  von  Willich  (Ref.  7),  carried  out  a 
Rayleigh-Ritz  solution  on  the  basis  of  a  variational  equation  equivalent 
to  the  governing  non-linear  differential  equations;  Chen  used  a  sixth 
degree  polynomial  (with  two  free  parameters)  for  the  representation  of  the 
deflection. 

The  results  of  the  more  recent  studies  of  Refs.  3-6  are  in  fairly 

good  agreement  in  the  narrow  range  between  X  =  U  and  X  =  5  j  beyond 

this  range,  the  results  do  not  encourage  confident  conclusions  concerning 

the  correct  theoretical  variation  of  p  with  X  . 

cr 

U.  FORMULATION  OF  INTEGRAL  EQUATIONS 

A  pair  of  integral  equations  equivalent  to  the  differential  equations 
(26)  and  (27)  together  with  the  boundary  conditions  (23)  and  ( 2U)  will  now 
be  constructed.  Consider  the  equations 


L(y)  +  xz  =  r 

(30) 

w 

II 

fr 

1 

tQ 

(31) 

wiiere  L(  )  is  the  operator  x(  ) 1  -  i(  )  .  Assume  that  the  functions 
r  and  s  vanish  for  x  >  X  ,  but  seek  a  solution  of  (30)  and  (31) 
valid  for  all  x  >  0  ,  and  bounded  at  infinity.  Application  of  the 
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Hankel  transform 


7  (a)  xJ1(ax)y(x)dx 


to  (30)  ani  (31),  and  use  of  the  relation  (Ref.  11) 


iSzir , .  a  * 


provides  the  transformed  equations 


2  x  k  /rN* 
-ay  +  z  =  (-) 


2  k  X  /  3  \K 
-a  z  -  y  =  (-) 


whence 


-a2&*  -  (£)* 
3fi  vXy 


*  (->* 
X  X 


Then  application  of  the  inversion  formula 


aJ^(ox)y  (a)da 


and  interchange  of  orders  of  integration  gives 


a.  \ 

y(x)  =  {  G(x,Or(4)d?  +  f  H(x,£)s(£)d£  (33) 


X  X 

z(x)  =  -  r  H(x,£)r(?)d£  +  f  G(x,£)sU)d£  (3U) 

o  o 


where  the  symmetrical  kernels  G  and  H  are 
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G(x,£) 


(-T — )  J-,(ax)J_  (a4.)d£ 
cT+1  1  1 


oo 

H(x,£)  =  -  f  (-J7— )  J,(ax)J  (at)df 
I  x  X 


Biese  infinite  integrals  are  evaluated  in  Appendix  A,  with  the  results, 

in  texms  of  Kelvin  functions, 

G(x,£)  =  bei'x  ker'f;  +  ber'x  kei'g  (0  =  x  =  £) 

<  <  (37) 

=  bei't;  ker'x  +  ber'i,  kei'x  (0  =  E,  =  x) 


H(x.£)  =  ber'x  ker'g  -  bei'x  W---' ' £  (0  =  x  =  £) 

<  <  (33) 

=  ber'g  ker'x  -  bei'g  kei'x  (0  -  %  =  x) 

Next,  with  a  view  to  tne  ultimate  enforcement  of  boundary  conditions 
at  x  =  X  ,  define  tht  functions  A^v*)  and  B,  (x)  as  the  solutions  of 
the  homogeneous  equations 

L(A)  +  xB  =  0  (39) 

L(B)  -  xA  =  0  (UO) 

and  the  non  homogeneous  boundary  conditions 

AjUS  =  o  (ia) 

XB^(X)  -  vB1(X)  =1  (h 2) 

Similarly,  let  A^(x)  and  B^(x)  be  the  solutions  of  (39)  and  (UO)  that 


satisfy 


^00  = 


XB^X)  -  vB2(X)  =  0 


Formulas  for  A^  ,  A^  ,  and  are  readily  found  in  terms  of  Kelvin 
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functions  (see  very  similar  calculations  in  Ref.  12)  and  are  displayed  in 
Appendix  B. 

The  desired  integral  equations  can  now  be  formulated;  it  can  be 
verified  that  a  solution  6  of  the  non-linear  system  (26),  (27),  (23), 
and  (2U)  satisfies  the  equations 

9(x)  =  y(x)  +  dA^(x)  +  eA2(x)  (U5) 

$(x)  =  z(x)  +  dB1(x)  +  eA^(x)  -  2px  (U6) 

where  y(x)  and  z(x)  are  given  by  (33)  and  (3U),  with 

r  =  0$  +  eh^  (1*7) 

s  =  -  |  02  -  eh6  (U8) 

and  where 

e  =  -y(x)  (1 *9) 

d  =  2pX(l-v)  -  [_Xz'(X)-vz(x)]  (50) 

This  formulation  otill  contains  a  single  differentiation,  in  Eq. 

(50)  for  the  scalar  d  .  But  note  that,  from  (31), 

(x2z')'  -  (xz)'  -  x2y  =  xs 

whence 


X 

X2z'(X)  -  Xz(X)  =  (x2y+xs)ax 
0 

Hence,  Eq,  (50)  may  be  replaced  by 

X 

d  =  2pX(l-v)  -  (l-v)z(X)  j  (x2y+xs)dx  (51) 

0 

The  integral  Eqs.  (U5)  and  (1*6)  constitute  the  basis  for  the  present 
analysis.  Two  different  viewpoints,  each  corresponding  to  well  defined 
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physical  conditions,  may  be  adopted  for  the  solution  of  these  equations. 
The  pressure  parameter  p  may  be  regarded  as  prescribed  and  monotonicaliy 
increasing,  and  solutions  for  0  and  §  sought;  this  specification  of 
the  pressure  corresponds  to  loading  the  shell  by  pressurized  air  in  a 
chamber  of  large  volume.  Alternately,  the  average  deflection  can  be 
specified,  and  a  solution  found  for  p  ,  as  well  as  0  and  $  ;  this 
approach  represents  pressure  loading  by  a  relatively  incompressible  fluid 
(say  mercury),  the  volume  of  which  is  under  control  at  all  times.  The 
average  deflection  is 


rw  dr 

from  which,  by  integration  by  parts  and  use  of  the  conditions 
w  =  dw/dr  =0  at  r  =  a  , 

w  =  f  x29  dx 
X  0 


The  convenient  average-deflection  parameter 


w  „  r,/-. 

p  =  ~2  H  =  2  L3(1"V  \ 


may  be  introduced;  then 


(52) 


(53) 


(5u) 


Figure  5  shows  hypothetical  curves  of  p  versus  p  for  the 
clamped  shallow  shell.  In  case  (a),  the  p  -  p  relationship  is 
one-to-one,  and  no  buckling  occurs;  in  cases  (b)  and  (c),  a  pressure 
monotonicaliy  increasing  from  zero  causes  buckling  as  soon  as  the  first 


local  maximum  of  the  curve,  p  ,  is  exceeded.  On  the  other  hand,  if  tne 
average  deflection  -were  assumed  to  increase  monotonically,  there  would  be 
no  sudden  snapping  of  the  shell  in  case  (b),  but  the  pressure  would  still 
reach  a  maximum  value  p  ,  and  then  decrease  smoothly.  In  case  (c ),  again 
under  prescribed,  increasing  average  deflection,  the  shell  would  tend  to 
snap  from  point  (l)  —  just  after  a  small  reduction  of  p  from  p  ^  —  to 
point  (?),  with  a  sudden  drop  in  pressure.  (The  designation  p  will  be 
retained  for  the  pressure  maximum,  even  though  snapping  does  not  occur  at 
Pcr  under  the  prescribed  volume  type  of  loading. )  The  most  precise 
determination  of  p  ^  in  cases  (b)  and  (c)  would  be  accomplished 
analytically  by  letting  p  increase  in  sufficiently  small  increments  and 
following  the  resulting  p  -  p  curve  beyond  its  first  maximum.  Alternatively, 
as  was  attempted  in  Ref.  J4.,  close  bounds  on  p^  can  bf.  deduced  by  letting 
p  increase  in  small  increments  until,  for  some  increment  in  p  ,  no 
"  adjacent"  equilibrium  position  is  found  to  exist.  Finally,  the  unlikely- 
looking  hypothetical  p  -  p  curve  shown  in  Fig.  5  cannot  be  ruled  out  on 
theoretical  grounds,  but  no  evidence  of  its  occurrence  was  found  in  any  of 
the  cases  actually  calculated.  (Note,  however,  that  such  curves  may  be 
found  when  the  pressure  is  plotted  against  the  deflection  at  some  particular 


point  in  the  shell  instead  of  against  the  average  deflection.) 


5.  ITERATIVE  PROCEDURE 

With  either  p  or  p  considered  prescribed,  an  iterative  procedure 
for  the  solution  of  Eqs.  (U5)  and  (U6),  can  now  be  formulated.  If  p  is 
prescribed,  an  initial  guess  for  8  and  may  be  substituted  into  the 
formulas  (33)  and  (3d)  for  y(x)  and  z(x)  ,  next  the  constants  e  and  d 
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can  be  calculated  by  Eqs.  (US')  and  (51) ,  and  then  the  subsequent 
approximation  to  0  and  ^  may  be  found  from  Eqs.  (U5)  and  (U6) .  The 
value  of  p  corresponding  to  the  new  Q  and  can  be  found  con¬ 
veniently  by  use  of  the  equation 


P 


"y  dx  +  dp  +  epr 


(55) 


where  the  constants 


(56) 


and 


\ 


x2A2 (x)dx 


(57) 


are  calculated  in  Appendix  B.  The  entire  process  may  then  be  repeated, 
and  a  point  on  the  p  -  p  curve  is  finally  determined  when  successive 
values  of  p  converge  to  sufficient  accuracy. 

If  p  is  prescribed,  y(x)  and  z(x)  are  found  in  the  same  way, 
and  so  is  e  ,  but  d  is  calculated  from  Eq.  (55),  and  Eq.  (51)  is  then 
used  to  solve  for  the  current  approximation  to  p  .  Again,  the  process 
is  repeated  until,  this  time,  the  values  of  p  converge  to  a  limit. 

In  the  application  of  these  procedures,  the  integrations  involved 
in  Eqs.  (33),  (3U),  (51)  and  (5U)  were  performed  numerically  on  the  basis 
of  Simpson's  rule,  and  the  iterative  procedure  was  formulated  in  terms 
of  matrix  operations.  (Special  integration  formulas  were  devised  for  the 
integrations  involving  the  kernel  G(x,£)  ,  which  has  discontinuous  first 


derivatives  at  x  =  1'  . ) 

The  spacing  of  stations  in  the  numerical  integration  is  very 
important;  in  this  work  the  invariant  spacing  ix  =  .25  was  chosen  for 
all  X  ,  leading  to  matrices  of  order  LiX  .  (Since  0  and  vanish 
at  x  =  0  ,  the  station  at  the  origin  could  be  omitted  in  the  matrix; 
set-up.)  This  choice  of  interval  was  motivated  primarily  by  con¬ 
sideration  of  the  waviness  of  the  expected  deformation  at  buckling;  the 
higher  X  ,  the  more  waves  are  expected,  with  a  roughly  constant  half- 
wave-length,  in  terms  of  x  ,  of  about  3  (see  Appendix  C).  The  use  of 
about  twelve  stations  per  half -wave -length  in  Simpson's  rule  was  con¬ 
sidered  adequate  for  acceptable  accuracy,  (it  happens,  also,  that  the 
choice  hx  =  .25  provides  a  reasonable  number  of  stations  within  the 
boundary  layer  of  early  bending  deformation. )  The  price  paid,  of  course, 
for  uniform  accuracy  over  all  values  of  X  was  the  necessity  of  handling 
matrices  of  increasingly  higher  order  as  X  increased.  The  use  of  high¬ 
speed  electronic  computation  was  essential;  an  IBM  70U  digital  computer 
was  used  to  obtain  the  numerical  results  presented  in  the  next  section.^ 

6.  NUMERICAL  RESULTS 

The  resuLts  obtained  for  the  variation  of  p  with  X  were 

cr 

established  by  a  combination  of  the  two  iterative  procedures  discussed. 

Fig.  6(a)  shows  sane  p  -  p  curves  for  the  initially  perfect  shell  obtained 
by  letting  p  be  the  prescribed  variable,  and  following  the  curve  over  the 
first  maximum  of  p  .  (Poisson's  ratio  v  was  taken  equal  to  1/3  in  all 
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The  author  is  indebted  to  Kr.  Morris  Waters  for  programming  the 
problem. 
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calculations.)  Points  on  these  curves  were  found  by  starting  with  a 
relatively  low  value  of  p  ,  tnen  letting  p  increase  in  increments  of 
.1,  and  finally,  when  necessary,  going  back  to  the  vicinity  of  with 

smaller  increments  of  p  (.01  -  « 02)  to  establish  p  more  accurately. 

In  these  calculations,  the  final  results  for  S  and  $  at  a  given  p 
were  used  as  the  initial  guesses  for  0  and  £  at  the  next  incremented 
value  of  p  .  At  each  p  ,  the  convergence  criterion  was  that  the  final 
answer  for  p  agree  to  three  decimal  places  with  the  mean  of  the  results 
of  the  five  previous  iterations. 

In  some  cases  (for  example,  1  =  5o5  ,  9  ,  e  =  0  ,  Fig.  6(b)),  the 
peak  of  the  p  -  p  curve  was  evidently  so  sharp,  and  the  subsequent 
decrease  in  p  as  well  as  p  so  sudden  (see  Fig.  5(c)),  that  increments 
of  .02  in  p  were  too  large  to  map  the  top  of  the  curve.  Thus,  for 
1  =  9,  the  p  -  p  variation  was  monotonic  up  to  p  -  1.3^  ,  but  at 
p  =  1.36  (shown  by  the  vertical  bar  in  Fig.  6(b))  the  iterative  pro¬ 
cedure  suddenly  failed  to  converge.  The  curve  was  then  corroborated  by 
use  of  the  other  procedure;  p  was  varied  in  increments  cf  .01,  and  at 
p  =  .9h  t  shown  by  the  horizontal  bar  in  Fig.  6(b),  no  convergence  for  p 
was  found.  The  loss  of  convergence  was  spectacular;  at  p  =  .93  *  p 
converged  after  12  iterations,  but  105  iterations  were  insufficient  at 
p  =  ,9U  .  The  value  of  pcr  for  X  =  9  ?  may  therefore  be  presumed  to 
be  less  than  .9U,  but  certainly  greater  than  .931,  the  value  obtained 
for  p  =  1.3h  . 

For  most  of  the  other  values  of  X  (and  e)  considered,  the  procedure 
based  on  assigning  p  was  omitted  altogether,  and  only  p  was  varied  in 
increments  of  .01  with  the  aim  of  establishing  sufficiently  cicse  bounds 
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on  pcr  „  In  most  cases,  the  sudden  failure  to  converge  to  a  definite  value 

of  p  at  some  value  of  p  was  as  dramatic  as  for  X  -  9  « 

The  results  found  for  p  at  all  the  values  of  X  and  e  considered 

cr 

are  given  in  Table  1„  Where  a  single  number  is  given  for  p^  s  the 
procedure  was  to  assign  p  and  pinpoint  accurately  the  top  of  the  p  -  p 
curve,.  Where  two  numbers  are  shown,  the  higher  value  corresponds  to  an 
unconverged  iteration  with  p  assigned,  and  the  lower  number  is  a  definite 
pre-buckling  value  of  p  obtained  by  prescribing  either  p  or  p  ,  as 
indicated  in  the  Table, 

The  results  are  shorn  in  Fig,  7  as  curves  of  p  versus  X  .  for 

cr 

various  values  of  e  „  The  crosses  on  these  curves  correspond  to  the  cases 
where  the  maximum  of  the  p  -  p  curves  were  precisely  located;  the  solid 
dots  show  the  points  for  which  only  bounds  on  were  found,  with  the 

diameter  of  the  dots  equal  tc  cr  greater  than  the  interval  of  uncertainty. 
Strictly  speaking,  it  must  be  conceded  that  the  upper  bounds  in  these  cases 
have  not  been  rigorously  established,  since  it  is  conceivable  that  failure 
to  converge  might  be  due  to  sane  unknown  cause  other  than  the  non-existence 
of  an  adjacent  equilibrium  position.  But  the  value  of  p  at  the  base  of 
each  dot  is  definitely  a  lower  bound  to  the  correct  value  of  p  , 

Also  shown  in  Fig,  7  are  experimental  points  obtained  by  Kaplan  and 
Fung  (Ref,  2)  and  Homewood,  Brine  and  Johnson  (Ref,  13),  The  significance 
of  the  present  theoretical  results  and  their  relation  to  these  experiments 
will  be  assessed  in  the  discussion  that  follows. 

7.  DISCUSSION  OF  RESULTS 

Comparison  of  Fig,  7  and  Fig,  U  reveals  that  only  for  X  5  5  are  the 
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present  results  for  the  initially  perfect  shell  (e  =  0)  in  fairly  good 
agreement  with  those  of  Archer*  Reiss  et  al,  Weinitschke,  and  Chen:;  for 
X  up  go  6,  Chen!s  energy  solution  remain/,  in  fair  agreement  with  the 
pieseni  one*  but  for  \  >  6  ,  the  present  results  bear  little  resemblance 
to  the  only  other  previously  available  ones  of  Archer  and  Reiss  et  al. 

In  common  with  the  results  of  Archer,  the  present  values  of  p  ^  appear 
tc  approach  1  as  X  increases,  but  do  so  in  an  oscillatoiy  fashion  rather 
than  monotonic ally  as  Archer  found. 

Some  explanation  for  the  details  of  the  oscillatory  variation  of 
Pci>  with  \  is  afforded  by  the  results  of  a  simple  analysis  of  the 
buckling  and  post-buckling  behavior  of  a  shallow  shell  having  boundary  con¬ 
ditions  that  are  slightly  different  from  those  of  complete  edge  fixity,. 

This  modified  shell  is  shown  by  the  upper  sketch  in  Fig,  8j  rotation  of 
the  edges  is  still  fully  constrained,  but  motion  normal  to  the  surface 
is  freely  permitted  at  the  boundary.  Thus,  just  as  the  complete  sphere 
does,  this  modified  shell  contracts  without  bending  under  uniform  pressure 
until  buckling  occurs.  The  scalloped  curve  in  Fig.  8  shows  the  results 
for  pcr  obtained  from  an  exact  solution  given  in  Appendix  Cj  each 
scallop  corresponds  to  the  lowest  ol  the  eigenvalues  in  the  spectrum  of 
critical  pressures  that  exists  for  overy  X  „  (Note  that  the  minimum  of 
each  scallop  is  at  p  =  1  . )  The  post-buckling  behavior  of  the  modified 
shell  is  also  analyzed  in  Appendix  C  by  means  of  a  simple  (but  approximate ) 
Galerkin  solutionj  typical  results,  for  four  values  of  X  ,  are  shown  in 
Fig.  9.  (At  any  X  ,  separate  post-buckling  curves  of  p  versus  p 
emanate  from  each  eigenvalue,  and  so  more  than  one  such  curve  can  often 
be  plotted  in  the  ranges  cf  p  and  p  shown.)  The  index  n  denotes 
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the  order  of  the  buckling  mode  shape,  and  is  numerically  equal  to  the 
number  of  nodal  circles.  The  shape  of  these  post-buckling  curves  indicates, 
of  course,  that  initial  imperfections  might  have  an  important  influence  on 
the  buckling  pressure  of  the  modified  shello  Notes,  however,  that  these 
curves  are  explicitly  based  on  the  assumption  that  the  center  buckling 
deflection  is  downward;  the  sign  of  the  eigenfunction  of  the  linear 
buckling  problem  is  indeterminate,  but  the  assumption  of  an  upward  buckle 
at  the  center  would  lead  to  p  -  p  curves  that  continue  to  rise  sharply 
after  buckling* 

The  shape  of  the  p^^  versus  X  curve  of  the  fully  clamped  shell., 
shown  for  comparison  in  Fig.  8,  can  now  be  rationalized.  Around  X  =  i;  , 
the  curve  is  roughly  similar  to  the  first  scallop  of  the  modified  shell, 
but  substantially  lower;  this  is  because  the  rotations  produced  by  the 
pressure,  having  the  distribution  shown  by  the  top  curve  in  Fig.  10(a), 
is  very  similar  to  the  buckling  rotations  (Fig.  10(b))  of  the  modified 
structure.  The  clamped  shell  is  predisposed  to  a  similar  buckling  mode, 
and  so  the  initial  rotations  produced  by  the  pressure  are  of  just  the 
right  kind  to  trigger  buckling.  Around  >  «=  7  ,  on  the  other  hand.  Fig.  10 
shows  that  the  early  rotations  could  hardly  be  less  sympathetic  to  the 
natural  buckling  rotations;  not°  that  the  sign  of  the  buckling  rotations, 
discussed  above,  is  important  in  this  comparison.  The  result  is  that 
buckling  is  delayed  until  p  reaches  values  somewhat  above  1,  the  initial 
distortions  being  of  little  or  no  help  in  inducing  buckling  in  the  pre¬ 
ferred  second  mode  shape.  (The  hesitant,  wavering,  approach  to  instability 
indicated  by  the  p  -  p  curve  for  X  =  7  (Fig.  6(a))  seems  to  reflect 
this  situation. )  A  happier  degree  of  compatibility  exists  between  the 
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third  mode  shape  and  the  early  bending  deformations  around  X  =  10  s  and 

in  this  region  the  third  mode  scallop  of  the  modified  shell  and  the 

Pcr  -  X  curve  of  the  clamped  shell  are  similar,,  Analogous  situations 

recur  as  X  increases,  but  already  at  X  =  13  the  initial  bending 

boundary  layer  (Fig.  10(a))  is  so  small  that  only  small  deviations  of 

p  from  1  are  found. 

*cr 

Not  unexpectedly,  the  curves  of  Fig.  7  corresponding  to  initially 
imperfect  shells  follow  the  same  trend  as  that  for  the  perfect  shell. 

(For  a  smaller  range  of  X  similar  results  —  not  directly  comparable, 
however,  because  of  different  choices  of  imperfection  shape  —  were 
calculated  by  Chen  in  Ref.  6  by  the  energy  method.)  The  reductions  in 
Pcr  caused  by  downward  imperfections  of  amplitudes  2  1/2  °/o  and  5  °/o 
of  the  shell  rise  are  substantial,  but,  except  near  X  =  k  ,  not  enough 
to  produce  even  fair  agreement  with  the  experimental  data  shown  in  Fig.  7. 
Furthermore,  imperfections  of  negative  sign  —  upward  with  respect  to  the 
perfect  shell  —  actually  increase  p^^  at  X  =  h  and  6. 

Kaplan  and  Fung  measured  the  initial  conto’irs  of  their  test 
specimens}  reduction  of  their  measured  imperfections  to  the  presently 
considered  form  (57)  reveals  that  they  had  values  of  e  that  averaged 
about  .15.  But  it  seems  evident  from  the  present  theoretical  results  that 
consideration  of  even  such  large  imperfections  would  still  not  fully  close 
the  gap  between  experiment  and  axisymmetric  buckling  theory  for  X  =  6  . 

The  implication  seems  inescapable  that  a  rational  explanation  of  the 
experimental  results  must  incorporate  consideration  of  non-axisymmetrie 
initial  imperfections  and  deformations.  The  modified  shell  has  a  spectrum 
of  eigenvalues  associated  with  buckling  modes  having  nodal  diameters  as  well 
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at>  circles;  for  X  greater  than  about  6,  the  clamped  shell  would  proi-ably 
choose  to  amplify  similar  initial  imperfections  ana  buckle  in  a  non- 
axisymmetrie  pattern.  In  oontrsist  to  the  situation  for  axisymmetric  im¬ 
perfections,  the  influence  of  an  imperfection  that  is  antisymmetric  through 
the  center  is  obviously  independent  of  sign.  It  can  only  be  speculated  now 
that  the  quantitative  influence  on  pc^  of  an  unsymmetrical  imperfection 
would  be  greater  than  that  of  an  axisymmetric  one  of  comparable  magnitude, 
but  this  seems  essential  if  the  experimental  results  are  to  be  explained. 
Finally,  there  is  experimental,  evidence  favoring  unsymmetrical  buckling  for 
X  >  6  .  Most  of  the  shells  tested  by  Homewood,  Brine  and  Johnson  snapped 
into  a  configuration  with  the  maximum  buckle  depth  off  center,  as  did  about 
half  of  the  specimens  of  Kaplan  and  Fung  in  the  range  1  >  6  .  Furthermore, 
even  a  final  buckled  shape  that  is  axisymmetric  might  be  reached  through  a 
succession  of  unsymmetrical  states,  as  has  been  shown  cle^arly  by  high  speed 
photographs  taken  by  Kloppel  and  Junpbluth  (Ref.  ill )  of  the  buckling  pro'  =ss 
of  some  non-shallow  spherical  shells. 


APPENDIX  A 


Calculations  of  Kernel  Functions  G  and  H 
Note  that  G(x,£)  as  given  by  Eq*  (35)  may  be  written 


G(x,0  =  (~r — -)Jq(ck)  (  Al) 

L0  a  +1  J 

Next.,  consider  the  Bessel  f  motion  of  the  second  kind  X^C2)  in  the  complex 
plane;  if  the  negative  imaginary  axis  is  taken  as  a  branch  line,  then  for 
x  real  and  positive 


V  { _  _  v  f-v-'N  _  0«.;T  /V'l 


(See,  for  example.  Ref.  15- )  Then,  the  integral  in  the  brackets  of  Eq.  (Al) 


is  given  by 


1  = "  n  (  (_r“)Jo(ax)  [Yi(a^  +  V-’0^] du 

0  a  +1 

00  2 

=-h  r  (-r“)jo(cx)Ti(^)dG 

J  a  +1 


where  the  Cauchy  principal  value  of  the  integral  is  understood.'  The  integral 


may  also  be  written 


op  2 

1  =  \  (  (^l)J0(ax)H{1)(a?)da 

J  au+l 


wherein  the  Hankol  f auction  of  the  first  kind 

=  Jx  +  iX1 


has  been  introduced. 


This  trick  ;ias  been  described  by  Sneddon  (Ref.  11,  p.  U76);  the  writer 
is  indebted  to  Professor  G.  F.  Carrier  for  suggesting  its  use. 
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The  asymptotic  behavior  of  the  integrand  (see  Ref.  15)  is  such  that 
for  0  -  x  -  £  ,  contour  integration  in  the  upper  half -plane  is  feasible 
for  the  evaluation  of  I  ,  and  the  result,  by  residue  theory,  is 


..1 


i2.JQ(xi2)Fj1)(€i2)  +  i  2J0(ad.3/2)H[1)(«i3/2)j 


(0  5  x  ^  K) 


(A3) 


In  terms  of  the  Kelvin  functions,  defined  by 

3/2 

ber  x  +  i  bei  x  =  JQ(xi  '  ) 
ker  x  +  i  kei  x  =  K^(xi^2) 

where  K  is  the  modified  Bessel  function  defined  by  K  (x)  =  -  ^  H^^ix) 
u  v  2  0 

the  result  (A3)  can  be  reduced  to 

I  =  bei  x  ker' ?  *  her  x  kei' 5  (0  =  x=  5) 

The  final  result  (37)  then  follows  by  differentiation,  and  the  symmetry 
of  G  . 

The  kernel  H  is  given  by 

BU,K)  -  § 


where 


oo 

L  =  f  (~r—  )J0(ox)J  (a^)da 
6  a+1 


(-J—)J0(ax)U^1\a^)da 

a+1 


L  =  ber  x  ker' 5  -  bei  x  kei’g  +  ~ 


Then 
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(tfcis  time  there  is  a  non-zero  contribution  from  the  semi-residue  at  the 
origin)  and  the  result  (38)  for  H  follows. 


*26- 


and  similarly  for  p, 
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&nd  similarly  (UU)  leads  to 


P2  ■  -  ^  (1-v)B2(\) 
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APPENDIX  G 

Analysis  of  Shell  with  Modified  Boundary  Conditions 
Foir  the  perfect ,  shallow  spherical  shell  supported  as  shown  by  the 
upper  sketch  in  Fig.  8,  the  membrane  state  6-0  and  £  =  -2px  is 
truly  an  exact  solution.  The  buckling  pressure  parameter  p  for  this 
modified  shell  is  found  by  setting 

$  =  -2px  +  0  (Cl) 

in  Eqs.  (19)  and  (20),  and  dropping  non-linear  terns  in  0  and  0  .  The 
resulting  linear  equations  are 

(x0')*  -  ^  +  x0  =  -2pxfi  (C2) 

(x0«)*  -  |  -  x9  -  0  (C3) 

The  boundary  conditions  are  0-0  and  0  =  0  at  x  -  X  ,  the  latter 
condition  following  from  the  equilibrium  requirement  §  (X)  -  -2pX  . 

The  spectrum  of  solutions  to  this  eigenvalue  problem  is 


where 

Ji(rn)  =0  (n  -  1,  2,  3,  ....)  (C6) 


and  the  corresponding  critical  pressures  are 


p 

Pn  -  |  (~)2  +  (x>2  <n  -  2,  3,  ....)  (C7) 

L  n 

The  roots  r&  of  (C5),  obtained  from  Ref.  16,  are  given  in  Table  2, 
below)  the  scalloped  curve  in  Fig*  8  shows  the  variation  of  the  lowest 
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critical  pressure  with  X  .  The  mode  shapes  (CU)  corresponding  to 
n  ■  1,  2,  3,  U  are  plotted  in  Fig.  10(b). 

The  post-buckling  behavior  of  the  modified  shell  may  be  deta rained 
approximately  by  means  of  a  simple  Galerkin  solution.  Substituting 


P  =  P  +  f 

(C8) 

0  - 

(C9) 

(CIO) 

r  x 

into  (19)  and  (20),  multiplying  these  equations  by  J^-S-) 

,  and 

integrating  from  zero  to  X  ,  gives 

,  9  ABXd 

(■“j  (A-B)  +  2Afn  + - =  0 

n  rn 

(Cll) 

td 

1 

> 

+ 

3°" 

n 

o 

(C12) 

where 

1  3 

2  f  [Vn^]  d« 
d"  *  [W]2 

(C13) 

The  constant  d  is  identical  with  Archer’s  constant  C1111 
n  n 


calculated  numerically  by  him  in  Ref.  3  ,  with  the  results  as  shown  in 
Table  2. 

Elimination  of  B  from  (GU)  and  (C12)  gives 

f  "  ^)2  ['3a  +  fl2]  (ClU) 

n 
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where 

“  -  (%)A 

The  average  deflection  is  given  by 

v  =>  w^  + 

where 

Pn  (C15) 

is  the  average  deflection  due  to  the  pure  contraction  of  the  shell  before 
buckling,  and  is  the  post-buckling  contribution  to  w  .  Then  (see 
Eq.  52) 

X 

W?  -  \v  ]  +  f  X20  dx  (C16) 

L  J  r=a  X4  ' 


Wl“ 


li(l-v  ;H 
L  X2  . 


But  note  that  the  modified  boundary  conditions  require  that 


[S*L  '  -  ®  M, 

where  u2  is  the  horizontal  displacement  during  buckling,  and,  by  (5), 

r  1  a  |>2  *2 

LU2j_  =  Et  Ldr  "  V  r. 

J  r=a  L  r*=a 


r»a 


where  if^  is  the  post-buckling  change  in  the  stress  function  f  .  Then 

.2 


2liH(l— v  ) 


(C17) 


where  ^2  =  §{x)  +  2p^x  is  the  change  in  $  after  buckling.  Then,  the 
non-dimensional  average-deflection  parameter  p  (Eq.  5L)  becomes,  from 
(CIS),  (06),  and  (C17) 


A. 

p  -  2(l-v)pn  -  [$2(X)  "  ^2^)J  +  "T  f  x20  dx 

X  o 

But  since  ©(x)  and  ^>(x)  should  satisfy  Eq„  (20),  a  manipulation  similar 
to  that  preceding  Eq.  (51)  of  the  text  may  be  applied,  giving 

X  X 

$2<X)  ■  £$2^  +  \  f  x20  dx  -  ~  ^  x©2  dx 
X  0  2X  n 

Vertical  equilibrium  now  imposes  the  condition  ^  (X)  =  -2pX  $  hence 
2(X)  -  -2X(p-Pn)  ,  and  so 


A. 

p  “  2(l-v)p  +  -i*  [  x02  dx 
■>\d  } 


(G18) 


Substituting  (CU)  into  the  integral  in  (Cl8)  gives 


p  -  2(l-v)p  +  l  [jQ(rn)]2  A2 


(C19) 


p-2(l-v)p  +  -^[j0(rn)]2a2 
an 

The  infomation  in  Eqs.  (C8),  (Clii),  and  (Cl?)  can  now  be  summarized 


as  follows: 


(a)  Before  buckling 


(b)  After  buckling 


(P  <  Pn  »  a  -  0) 
p  =  2(l-v)p 


P  *  Pn  ♦  r  ( p-)2(-3 a+a^) 
n 

p  -  2(l-v)p  +  J  (£-)2  bn  a2 
n 


(C21) 
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where 


b 


n 


2 


r 

n 


7 


2 


n 

The  constants  b  are  tabulated  in  Table  2* 
n 

The  post-buckling  variation  of  p  with  p  is  given  parametrically 
by  Eqs.  (C21).  Typical  p  -  p  curves  (with  v  =  i  )  are  given  in  Fig.  9 
for  several  X's  ,  far  a  >  0  .  Note  that  a  <  0  would  give  quite  different 
curves,  involving  no  reduction  in  p  j  hence  it  may  be  presumed  that  a  >  0 
—  meaning  downward  buckling  at  the  center  —  is  always  preferred. 
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TABLE  1 

BUCKLING  PRESSURES 


(a)  Perfect  Shell  (b)  Imperfect  Shell 


X 

D 

Prescribed 

Prescribed 

*cr 

Parameter 

£ 

Pcr 

Parameter 

3.5 

.6lU 

P 

1* 

—  025 

.617 

P 

h 

.578 

P 

1* 

-.050 

.659 

P 

5 

.62  9 

P 

1* 

.025 

1.061 

P 

5.5 

.789-. 19 

P»P 

k 

.o5o 

1.116 

P 

6 

.: 99 5 

P 

6 

.025 

.90-. 91 

P 

7 

1*068 

P 

6 

.o5o 

.62-. 63 

P 

8 

1.13C 

P 

8 

.025 

1.03-1.0U 

P 

9 

.931-. 9h 

P»P 

8 

.o5o 

.93-91* 

P 

10 

.822-8 3 

P»P 

10 

.025 

.78-79 

P~ 

11 

.83-81* 

P 

12 

.025 

.83- 8k 

P 

12 

.96—97 

P 

12 

.o5o 

.75-76 

P 

13 

.96-. 97 

P 

TABLE  2 

CONSTANTS  IN  ANALYSIS  OF  MODIFIED  SHELL 


n 

r 

n 

d 

n 

b 

n 

1 

3.83 

1.012 

2.33 

2 

7.02 

0 

CO 

H 

CO 

6.62 

3 

10.17 

.901 

7.95 

1* 

13.32 

.81*6 

11.82 

5 

16.1*7 

.883 

13. kl* 

Note;  For  large  n  ,  d^^/2  ,  bn -*  (8/3n)rn 


FIG.  1  CLAMPED  SHALLOW  SHELL  OF  REVOLUTION 


PRESSURE 


FIG.  2  PRESSURE -DEFLECTION  CURVES  FOR 
PERFECT  (A)  AND  IMPERFECT  (B) 

(C)  COMPLETE  SPHERICAL  SHELLS 


FIG.  3  DISPLACEMENTS,  ROTATION,  FORCES  AND  MOMENTS 
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FIG.  4  PREVIOUS  RESULTS  FOR  BUCKLING  PRESSURES 
OF  CLAMPED  SHALLOW  SPHERICAL  SHELLS 


p  p 
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FIG.  5  HYPOTHETICAL  CURVES  OF  PRESSURE 
VS.  AVERAGE  DEFLECTION 


FIG.  7  CALCULATED  BUCKLING  PRESSURES  OF  INITIALLY  PERFECT  AND 
IMPERFECT  CLAMPED  SHALLOW  SPHERICAL  SHELLS  AND 
EXPERIMENTAL  RESULTS 


FIG.  8  COMPARISON  OF  AXISYMMETRIC 
CLAMPED  SHALLOW  SPHERICAL 
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FIG.  10  COMPARISON  OF  (a)  EARLY  (LINEAR  THEORY)  ROTATIONS  OF  CLAMPED  SHELL 
AND  (b)  BUCKLING  ROTATIONS  OF  MODIFIED  SHELL 
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